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In view of new constraints put forth by recent observations and measurements in the realm 
of astrophysics and nuclear physics, we update the non-linear realization of the sigma model as to 
reflect such constraints. By doing this, we obtain new equations of state that may be used to describe 
neutron stars. Such equations of state are obtained by investigating different ways by which the 
vector mesons self-interact. Furthermore, we also investigate the role played by the delta mesons in 
the model. As a result, we are able to develop equations of state that are in better agreement with 
data, such as nuclear compressibility and slope of the symmetry energy at saturation, star masses, 
radii, and cooling profiles. 


In recent years, measurements of massive neutron stars 
with masses around 2.0 M© mm have been used to select 
equations of state that fulfill astrophysics constraints. In 
addition, estimates of neutron star radii, which seem to 
be smaller than previously thought, have become much 
more sophisticated mm- These have been used with the 
same purpose. Note that, in addition to more direct radii 
measurements, trustworthy theoretical calculations can 
provide insight into this quantity [MI]. Still, until more 
accurate radii measurements are made, it is important to 
explore a somewhat large range of values. 

Nuclear physics constraints have also improved, and 
quantities such as the symmetry energy and its slope have 
had their values narrowed at saturation density to 29.0 
MeV < Tisym < 32.7 MeV and 40.5 MeV < L < 61.9 
MeV in the analysis of Ref. [12], which we take as a 
guideline for the underlying investigation. Other works 
have also found similar values |13j . For a recent review 
on the symmetry energy, please see Ref. M- 

In view of the aforementioned improvements in the 
analysis of astrophysical and nuclear data, specifically the 
slope of the symmetry energy at saturation density, we 
present an updated version of the hadronic chiral SU(3) 
sigma model. Following the idea proposed recently in 
Refs. [TSHTT] . we study different possible self-interactions 
of the vector mesons. We do this by considering differ¬ 
ent possible coupling schemes, but restricting ourselves 
to only those that are chirally invariant. 

We also investigate the influence of the delta meson 
on the model and on its predicted microscopic quantities 
and macroscopic observables. We recall that effects of 
the delta meson in relativistic mean held models have 
been investigated in Refs. and many others. 

Chiral sigma models are effective relativistic models 
that describe hadrons interacting via meson exchange 
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and, most importantly, are constructed from symmetry 
relations. They are constructed in a chirally invariant 
manner since the particle masses originate from inter¬ 
actions with the medium and, therefore, go to zero at 
high density and/or temperature. The underlying ex¬ 
tended sigma model is in very good agreement with nu¬ 
clear physics data, such as the vacuum masses of the 
baryons, saturation properties, hyperon potentials, and 
pion and kaon decay constants and A, etc. [2^1^ . 

The Lagrangian density of the SU(3) non-linear sigma 
model, called hereafter the chiral mean-field (CMF) 
model, applied to neutron star matter can be found in 
Refs. [211 di]. A recent extension of this model also in¬ 
cludes quarks as dynamical degrees of freedom [2SHS0]. 
In this work, we make use of a simple hadronic version of 
the model without hyperons or quarks. The reason be¬ 
hind such a choice is that we want to study the effect of 
different possible self-interactions of the vector mesons, 
without the interference of the appearance of hyperons 
or quarks in the star. The CMF Lagrangian density of 
the model we use in this work reads 


-C = £Kin + >Clnt + £self + >CsB , (1) 

where, besides the kinetic energy term for hadrons and 
leptons (included to ensure charge neutrality), the terms: 


Cint = - X V'*[7o(5ituW -b -b g^pT3p) + M* ]-0i, 

i 

^Seif = + mlp^ -b 




— k/^ 


In 


o'oCo 


“b Lvec4 


5^)C 


C. 


SB 


= - (^V2mlfk - 


C, (2) 




2 


represent the interactions between baryons and vector 
and scalar mesons, the self interactions of scalar and vec¬ 
tor mesons, and an explicit chiral symmetry breaking 
term, which is responsible for producing the masses of the 
pseudo-scalar mesons. The index i denotes the protons 
and neutrons. The electrons are included as a free Fermi 
gas. The mesons included are the vector-isoscalars u and 
(j) (strange quark-antiquark state), the vector-isovector 
p, the scalar-isoscalars a and (strange quark-antiquark 
state) and the scalar-isovector 6, with ra being twice the 
isospin projection of each particle (±1). The isovector 
mesons affect isospin-asymmetric matter and thus, are 
important for neutron star physics. Also, the 6 meson 
has a contrary but complementary role to the p meson, 
much like the a and w mesons. 

Note that the strange mesons influence the model 
through the coupling to the nucleons (C) and through 
self-interaction terms {(p). The different possible self¬ 
interaction terms of the vector mesons that are chiral 
invariant (referred to as Tvec 4 above) are the following 

• /:vec4 = 2Tr(y4) , 

. /:vec4 = [Tr(l/2)]2 ^ 

. /:vec4 = [TriV)]y4, 

where V stands for the matrix of the vector meson mul- 
tiplet, which reduces to a diagonal form in the mean field 
limit, i. e., 

/{uj + p)/V2 0 0\ 

v=\ 0 {UJ- p)/V 2 0 . (3) 

Vo 0 pj 

The different self-interaction terms of the vector mesons 
shown above correspond, respectively, to the coupling 
schemes Cl, C3, and C4 in the following. We also in¬ 
troduce the coupling scheme C2, which is a linear combi¬ 
nation of 1 and 3, and presents no ujp mixing. Note that 
the coupling scheme C4 is quite special, as it contains 
contributions that are linear in the isoscalar vector field 

• Cl: £vec4 = ff4(w"^ -I- Guj'^p'^ + p"^ + 2(f)'^) , 

• C2: £vec4 = ff4(w"^ + p'^ + /2 -I- ip^cp? -I- iuP'pp-) , 

• C3: £vec4 = 94 ( 1 ^^ + 2uj^p^ + p'^ + + 

2p^P^) , 

• C4: £vec4 = <74(w^ + + 

The effective masses for the baryons are simply gen¬ 
erated by the scalar mesons, except for a small explicit 
mass term Mq 

+ 5<C + ^ 0 . • (4) 

The coupling constants of the model were modified 
from Ref. [23], in the sense that we used ranges, rather 
than values, for nuclear constraints at saturation, such 



FIG. 1. (Color online) Slope of the symmetry energy at sat¬ 
uration as a function of the coupling constant of the <5 meson 
for different coupling schemes. 


as baryon density (po = 0.15 to 0.17 fm“^), binding en¬ 
ergy {BjA = —15.75 to —16.3 MeV), and symmetry en¬ 
ergy {Egym = 29 — 32.7 MeV). Within this setup, we try 
to achieve the minimum possible compressibility K and 
symmetry energy slope L, while still fulfilling astrophys- 
ical constraints. For a review on compressibility values 
see Ref. [3T] and references therein. 

More specifically, we vary the vector sector of the 
model, 94 , and Pp {g^ = 0) while keeping most of 
the the scalar sector constant {g^ = —9.83, = 1.22, 

fco = 1.19x^, fci = —1.40, ^2 = 5.55, fca = 2.65% and 
k 4 = —O.OOy"', with y = 401.93 MeV). The scalar sector 
of the model cannot be further adjusted as it was fixed to 
reproduce the vacuum masses of the baryons and scalar 
mesons, and the pion and kaon decay constants and 
/«. Only gs is modified freely, as it only affects isospin 
asymmetric matter. 

We now present microscopic and macroscopic proper¬ 
ties reproduced by the 4 coupling schemes Cl to C4. We 
fit the coupling constants 94 and g^i to reproduce nuclear 
and astrophysical constraints. We set gp so as to obtain 
two curves, one for the lowest accepted value from the 
range of the symmetry energy (29 MeV) and, another 
for the corresponding highest value (32.7 MeV). Finally, 
we study the behavior of different quantities as a func¬ 
tion of the strength of the delta meson coupling to the 
nucleons gs- 

Within this setup, it is interesting to see the behavior 
of L (the slope of the symmetry energy at saturation) as 
we vary the delta coupling. As shown in Fig. although 
the value of L increases with the delta coupling for cases 
C2, C3, and C4, it reaches a minimum value for case Cl 
around gs = 11.5. This is related to the strength of the 
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FIG. 2. (Color online) Maximum possible star mass as a 
fnnction of the coupling constant of the S meson for different 
coupling schemes. 


top term in the vector-meson self-interactions (present in 
Cl and C3 cases), which induces a softer equation of 
state for asymmetric matter at higher densities (due to 
lower absolute values of uj and p necessary to fit nuclear 
saturation and asymmetry properties). Note that, even 
when gs = 11.50, the coupling strength of the delta me¬ 
son is still small compared with that of the sigma meson 
{gslms)’^l{galmcrY ~ 0-3- As expected, L is higher for 
higher values of the symmetry energy for all coupling 
schemes, as pointed out in Refs. [32H35] . where a linear 
behavior between L and iilsym was found. 

The respective maximum masses achieved for stars 
generated by each coupling scheme are calculated 
through the use of the Tolman-Oppenheimer-Volkoff 
(TOV) equations and are shown in Fig.[^ Once more, the 
strength of the Lop coupling plays an important role, as it 
renders the equation of state softer, which in turn yields 
less massive stars. Nevertheless, the maximum masses 
for all the coupling schemes analyzed increase with gg^ 
since this isovector channel generates a larger asymmetry 
between protons and neutrons (at low and intermediate 
densities) and, therefore, more Fermi pressure in the star. 

As already mentioned, the coupling scheme C4 for the 
self-interaction of the vector mesons does not follow the 
trend set by the other schemes, in the sense that it con¬ 
tains a term which is linear on the strange vector meson 
(j). This feature modifies the behavior of the model, which 
has to be treated separately from the previous coupling 
schemes. For example, it has a bare mass for the nucleons 
of Mq = 50 MeV. This allows for smaller compressibil¬ 
ity (which otherwise would be large) in better agreement 
with nuclear physics data. Note that this version of the 
vector-meson self-interaction fitted to only one value of 



FIG. 3. (Golor online) Radius of the maximum possible star 
mass as a function of the coupling constant of the S meson 
for different coupling schemes. 


symmetry energy has already been discussed in Refs. Ei- 

EHES]. 

In order to discuss the radii of stars, it is important 
to include a separate equation of state for the crust. In 
this work, we choose to use the Baym-Pethick-Sutherland 
(BPS) equation of state containing an inner crust, an 
outer crust and an atmosphere [55]. The respective radii 
for each maximum mass star (shown in Fig. for each 
coupling scheme) are shown in Fig. It is clear that 
the radii behave in the opposite manner with respect to 
the increase of gs'. the radius increases and is in worse 
agreement with data [3H6]. Note that, for the coupling 
scheme C4, neither the maximum mass nor the respec¬ 
tive star radius increases much with gs, when compared 
with the other coupling schemes. This is again related to 
change of the value of the meson fields (specifically the 
lack of change in the meson cr, the counterpart of lo) with 
gs, which affects the equation of state. If, instead of radii 
for maximum mass stars we had shown the radii for 1.4 
Mq stars, Fig. [^ would look similar, except that all radii 
would be shifted up by about 2 km. 

At this point, it becomes clear that the influence of the 
S meson on isospin asymmetric equations of state has to 
be carefully analyzed. This becomes even more evident 
when we look at Fig. [^ which shows the increase of a 
star’s possible maximum mass and respective radius in 
relation to the different coupling schemes. At first, it 
looks like there is no reconciliation between large mass 
and small radius stars. However, after careful analysis, 
we can identify a parameter set that fulfills all the con¬ 
straints discussed. Such a parameter set is shown in Ta- 
ble[T| Note, in particular, that the parametrization chosen 
for Cl reproduces small slopes for the symmetry energies 
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FIG. 4. (Color online) Maximum possible star masses and 
respective radii for each coupling schemes obtained by varying 
gs from 0 to 15. 


(for the whole range of symmetry energies at saturation) 
and yet still yields massive stars. 

Finally, as a complementary approach to the usual 
study of neutron star macroscopic properties, we briefly 
discuss their cooling. The cooling or thermal evolution 
of neutron stars is strongly dependent on whether or not 
the direct Urea process (DU henceforth) takes place, as 
well as on pairing among the nucleons. Given efficient 
neutrino generation in the DU process, if it takes place 
in the star, it will lead to a fast (sometimes also called 
enhanced) thermal evolution of the object. This in turn, 
may lead to disagreement with observed thermal data 
on compact objects. The efficiency of the DU process 
highlights the need for an accurate description of pair¬ 
ing among nucleons in neutron stars, since it suppresses 
the neutrino emission from the DU process (among other 
processes as well) and may thus reconcile theory with ob¬ 
servation. For a detailed review on the neutrino emission 
of neutron stars as well as of pairing see Ref. [37] . 

Figure shows the density at which the DU process 
conditions are achieved, thus signaling its onset. This 
density is much higher (and the process consequently 
only takes place much more towards the center of the 
star) for the coupling schemes Cl and C3. These are 
again the cases in which there is a strong ujp term for 
the vector-meson self-interaction. Also, the DU thresh¬ 
old density decreases quickly with the increase of gs, as 
the star becomes more isospin symmetric (more protons 
at large densities) due to lower absolute values of the p 
meson. Finally, note that larger symmetry energies are 
associated with lower density thresholds. This is only 
natural, since the symmetry energy is directly related to 
the energy balance between symmetric and neutron mat¬ 



FIG. 5. (Golor online) Normalized density threshold for the 
direct Urea process as a function of the coupling constant of 
the S meson for different coupling schemes. 


ter. In other words, larger symmetry energy means larger 
energy cost to create more neutrons and, consequently 
larger proton fraction. The relation between symmetry 
energy and neutron star cooling has been studied, for 
example, in Ref. |55fHT] . 

To complement the previous discussion, we perform a 
preliminary study of the thermal evolution of neutron 
stars using the coupling schemes presented in this pa¬ 
per. More precisely, we performed full thermal evolution 
simulations of neutrons stars whose inner compositions 
are described by the parameter sets outlined in Table 
All relevant thermal process are consistently taken into 
account. At first, we do not consider any pairing in or¬ 
der to be able to clearly identify the effects of each cou¬ 
pling scheme on the thermal evolution. For this case, we 
show in Figs. and the thermal evolution of stars with 
M = 1.4 Mq and the maximum mass of each respective 
parameter set. 

Figures and show the effects of the onset of the 
DU process on the thermal evolution of the stars. As 
can be seen in Fig. the coupling scheme that exhibits 
the “slowest” cooling is C3 with a symmetry energy of 29 
MeV. This stems from the fact that for stars with that 
mass the DU process onset is not achieved within this 
parameter set (in agreement with Fig. [^. As for stars 
of the maximum mass, the DU onset is achieved in every 
model, as shown in Fig. 

Finally, to illustrate the importance of pairing, we per¬ 
form cooling simulations allowing the nucleons to form 
pairs. The neutrons are allowed to form singlet pairs 
in the crust and triplets in the core, whereas the pro¬ 
tons are allowed to form singlet pairs in the core. A 
similar description has been used to analyze the thermal 
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TABLE I. Best fit for the coupling constants of the model using different vector-meson self-interactions which fulfill nuclear 
physics and astrophysics constraints. The radius shown is the one corresponding to the star’s maximum mass. 


c 

54 

Quj 

9p 

as 

Mo (MeV) 

K (MeV) 

(MeV) 

L (MeV) 

M (Me) 

R (km) 

Cl 

58.40 

13.66 

11.06 

11.50 

0 

271.41 

29.00 

54.00 

2.05 

12.56 

Cl 

58.40 

13.66 

11.33 

11.50 

0 

271.41 

32.70 

59.05 

2.04 

12.49 

C2 

58.40 

13.66 

3.51 

0 

0 

271.41 

29.00 

85.68 

2.00 

12.22 

C2 

58.40 

13.66 

4.02 

0 

0 

271.41 

32.70 

96.77 

2.01 

12.32 

(73 

58.40 

13.66 

3.82 

0 

0 

271.41 

29.00 

77.19 

2.12 

14.38 

(73 

58.40 

13.66 

4.38 

0 

0 

271.41 

32.70 

85.64 

1.95 

12.01 

(74 

34.20 

11.88 

3.70 

0 

150 

292.91 

29.00 

83.38 

2.14 

12.18 

(74 

34.20 

11.88 

4.19 

0 

150 

292.91 

32.70 

94.48 

2.15 

12.26 



FIG. 6. (Color online) Surface temperature as observed at 
infinity of neutron stars with 1.4 M© for the parameters de¬ 
scribed in Table m 



Age (years) 

FIG. 7. (Color online) Same as Fig. but for neutron stars 
with the maximum possible mass. 


evolution of the neutron star in the supernova remnant 
Cassiopeia A [3^|33]. We note that, in addition to pair¬ 
ing, we have also included pair breaking-formation (PBF) 
neutrino emission process. 

Fig. [8] shows the thermal evolution of neutron stars in 
which pairing is present, as discussed above. In this case, 
the cooling of neutron stars with 1.4 Mq is very similar 
for all parameter sets and falls within the dark band in 
the figure. For neutron stars of higher masses we have 
two different scenarios, a “faster” cooling for the higher 
symmetry energy parameter sets of the coupling schemes 
Cl and C2 and both parameter sets of C4, and a “slower” 
cooling for the lower symmetry energy parameter sets of 
the coupling schemes Cl and C2 as well as both parame¬ 
ter sets of C3. These results are once more in agreement 
with the results shown in Fig.j^if we simply consider that 
the thresholds for the DU process are effectively raised 
by the presence of pairing. 

In conclusion, the aim of this work was to update the 
non-linear realization of the hadronic sigma model to 
reflect new astrophysical observations and experimental 
nuclear data. More specifically, we looked for equations 


of state that reproduce small star radii (at the same time 
as large star masses) and a small symmetry energy slope 
at saturation. In addition, we verified that the thermal 
evolution of stars reproduced by our equations of state 
were in agreement with Cas A cooling data. The main 
features of these equations of state are provided in Table 
I and all of our results (in the format of data tables) can 
be provided upon request. 

More specifically, we investigated the role of the vector- 
meson self-interaction terms within a chiral mean field 
formulation by performing a parameter scan. We also 
varied the coupling of the scalar-isovector meson to the 
nucleons. We found that the strength of the ujp term in 
the vector-meson self-interactions plays an important role 
in the equation of state, together with strange-vector- 
meson linear terms. 

With respect to neutron star thermal evolution, our 
preliminary study showed that our results are in good 
agreement with modern cooling calculations. In partic¬ 
ular, the schemes that include strong vector-meson self¬ 
interaction terms with low scalar-isovector coupling and 
the parametrizations that reproduce lower symmetry en¬ 
ergies have an advantage over the others, as they exhibit 
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FIG. 8. (Color online) Combination of Figs andwith the 
inclusion of pairing effects. The thermal evolution of 1.4 Mq 
stars from all parameter sets falls within the dark band. The 
thermal evolutions of the maximum mass stars are indicated 
by individual curves. 


a “slower” cooling that can be identified for lower mass 
stars. This effect is further enhanced when pairing is in¬ 
cluded in the calculations, which was done by following 
the prescription used to model the thermal evolution of 
Cas A. In this case, the difference between “slow” and 
“fast” cooling stars was significantly enhanced. A fur¬ 
ther, more in depth study of the effects of pairing in the 
thermal evolution of neutron stars using our equations of 
state will be reported in a future work. 
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